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The status of lattice determinations of quark masses is reviewed (with the exception of mb). Attempts to 
extract the low-energy constants in the effective chiral Lagrangian are discussed, with special emphasis on those 
couplings which are required to test the hypothesis of a massless up-quark. Furthermore, the issue of quenched 
chiral logarithms is addressed. 



1. INTRODUCTION 

Recently, many studies of QCD at low ener- 
gies have investigated the interplay between lat- 
tice simulations and Chiral Perturbation Theory 
(ChPT). I shall start the review of these activities 
by recalling the basic features of ChPT and point 
out what can be learned from combining ChPT 
with lattice simulations. 

Chiral Perturbation Theory is a systematic ex- 
pansion in the 4-momentum p and the quark 
masses m u , md, m s about the massless limit [Q. 
In the low-energy regime its information content 
is equivalent to that of QCD, and this has been 
exploited in many applications, e.g. in the cal- 
culation of pion scattering amplitudes and quark 
mass ratios. Furthermore, ChPT provides valu- 
able input for lattice simulations. Here, the pre- 
diction of the quark mass dependence of rrip S 
is surely the most widely used piece of informa- 
tion. In addition, ChPT can model the volume 
dependence of observables, and also relates dif- 
ferent processes. For instance, the amplitudes of 
K — » 7T7T decays can be expressed in terms of the 
theoretically much simpler K — * tt transition. 

However, Chiral Perturbation Theory is an ef- 
fective, non-renormalisable theory, parametcrised 
in terms of a set of empirical couplings, which are 
usually called "low-energy constants" (LECs). At 
lowest order the effective chiral Lagrangian reads 

+B Jr (M(U+tf))}, (1) 

where 

U = exp{iT a < / P/(y /2F )} (2) 

* Plenary talk presented at Lattice 2002, MIT, Cambridge, 
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is the field of Goldstone bosons, and M. = 
diag(m u , rrid, m. s ) is the quark mass matrix. The 
LECs at leading order are Bq and Fq, where the 
latter corresponds to the pion decay constant in 
the chiral limit ^ At next- to-leading order there 
are 12 more interaction terms and hence 12 new 
LECs, L\, . . . , L\2- 

The values of the LECs can be fixed by using 
experimental data and comparing with the rele- 
vant expressions obtained in ChPT. It turns out, 
however, that the complete set of LECs cannot be 
determined in this way. In other words, the val- 
ues of some of the LECs cannot be fixed without 
resorting to additional theoretical assumptions. 
One particular example is the value of B , which 
appears in the chiral expansion of the pion mass 
at leading order: 

m w = 2B m, m — \ {m n + md). (3) 

This shows that Bq can only be determined from 
m w if the physical value of the quark mass fh is 
known already. By the same token, the quark 
mass fh can only be inferred if an estimate for 
Bq is available. However, Bq drops out in suit- 
ably chosen ratios of mj, m^, .... Thus, while 
ChPT enables us to compute quark mass ratios, it 
fails to provide an absolute normalisation of their 
masses. Another reason why the complete set of 
LECs cannot be determined from chiral symme- 
try considerations alone is the fact that the ef- 
fective chiral Lagrangian beyond leading order is 
invariant under a symmetry transformation in- 
volving the LECs and the mass matrix Ai . This 
is the famous "Kaplan-Manohar ambiguity" [0 . 
At this point it is clear that lattice simulations 

2 Throughout this review I use conventions where F n = 
93.3 MeV. 
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can in turn provide valuable input for ChPT. By 
studying the quark mass dependence of Gold- 
stone bosons in a simulation one can determine 
the value of B , or, equivalently provide an ab- 
solute normalisation of quark masses. Since the 
Kaplan- Manohar (KM) ambiguity is not a sym- 
metry of QCD, it is possible to resolve it by de- 
termining the values of LECs from an approach 
based on first principles. It must be kept in mind, 
though, that a successful combination of lattice 
QCD and ChPT requires sufficient overlap be- 
tween the range of quark masses used in simula- 
tions and the region of validity of the chiral ex- 
pansion. This requirement is crucial if ChPT is 
used to extrapolate results from simulations per- 
formed for relatively heavy quarks and the regime 
of the physical u- and d-quark masses || . 

The remainder of this review covers a sum- 
mary of the current status of lattice determina- 
tions of the light quark masses, quenched chiral 
logarithms, and LECs at NLO. 

2. QUARK MASSES 

Recent years have seen a great deal of progress 
in lattice determinations of quark masses. Some 
of the dominant systematic effects are now un- 
der much better control, chiefly due to using 
improved lattice actions and non-perturbative 
renormalisation, as implemented cither via the 
RI/MOM scheme Q or the Schrodinger func- 
tional (SF) ||. Several collaborations have 
started to quantify quenching errors, by perform- 
ing simulations with dynamical quarks. Finally, 
the results in the quenched approximation are 
also being tested in simulations employing exact 
chiral symmetry on the lattice. Here we will focus 
on determinations of m s , fh and the charm quark 
mass. The mass of the &-quark is discussed in 
R. Sommer's contribution to this conference ||. 
Earlier reviews can be found in refs. 

2.1. Wilson and staggered fermions 

The LEC Bq is obtained straightforwardly by 
mapping out the quark mass dependence of pseu- 
doscalar meson masses. In order to compute, say, 
m s + fh one then has to specify the lattice scale 
and the quantity which fixes the bare physical 
quark mass. For instance, using ro to set the scale 



[l(J and the kaon mass m-K ( "if -input" ) one has 



i(m s + m)r = Z- 



B r 



x (™KTo)exp-( 4 ) 



Here, Z denotes the renormalisation factor which 
relates the quark mass in lattice regularisation 
to a reference continuum scheme (e.g. MS), and 
(fiKTo)exp * s t ne phenomenological estimate of 
the kaon mass in units of the lattice scale. In 
order to estimate m s , one either has to extract fh 
separately through an extrapolation to (m w ro) 2 , 
or take the ratio m s /fh from ChPT. 

A compilation of recent results for m s and the 
ratio m s /fh is shown in Table |l]. A marked feature 
is that estimates for m s in the continuum limit 
have stabilised, in contrast to the situation seen 
a few years ago. The table also shows that lattice 
estimates for the ratio m s /fh, which are all based 
on chiral extrapolations, are broadly consistent 
with the result in ChPT at NLO 1531: 



m s /fh = 24.4 ± 1.4. 



(5) 



However, closer inspection reveals small but sig- 
nificant deviations among the results for m s . 
Therefore, for the results to be consistent, these 
deviations must be due to quenching effects, 
which are caused by different choices for the lat- 
tice scale. We will now attempt to convert some 
of the results in Table [l] to a common scale, in 
order to check whether consistency is satisfied. 

To this end we consider two lattice scales, Q 
and Q' . From the definition of the quark mass 
in eq. ^ it follows that the strange quark mass 
mi in MeV, estimated using Q, is related to its 
counterpart mi® via 



,(Q) 



[MeV] 



Q' 



Here, the subscripts "lat" and "exp" refer to lat- 
tice and experimental estimates of the scale ra- 
tio, respectively. The deviation of the factor F 
from unity is indicative of the relative quench- 
ing effects, when either Q or Q' is chosen to set 
the scale. I have taken results for vector meson 
masses from refs. [[tl 12,17| and for Fk from 
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Table 1 

Results for m s in the MS-scheme at fj, = 2 GeV and for the ratio m s /m, obtained using Wilson (W) and 
staggered (KS) quarks (K- input). Where applicable, we include information on the implementation of 
0(a) improvement (non-perturbative (NP), mean- field perturbative (MF), tree-level (tree)), as well as 
the method for non-perturbative renormalisation (RI or SF). 



Collaboration 



Ref. 



[12| 
[13| 
[14| 
[15| 
[16| 
[17| 
[18| 
[19| 



action Impr. a [fm] Scale Ren. m s [MeV] m s /fh 



SPQ cd R 

CP-PACS 

CP-PACS 

APE 

QCDSF 

ALPHA /UKQCD 

JLQCD 

APE 

GGRT 



W 
W 
W/Iwas. 
W 
W 
W 
KS 
W 

w 



NP 

MF 
NP 
NP 
NP 

NP 

tree 







0.07 




0.07 
0.07 



771 if. 

m p 
m p 

ro 
F K 
m p 

m K * 
m K * 



RI 



RI 
SF 
SF 
RI 
RI 
RI 



106(2)(8) 
114(2)(1) 

no(l) 

111(9) 
105(4) 
97(4) 
106(7) 
111(12) 
130(2)(18) 



24.3(2)(6) 

26.5(|;1) 

25.0(1.4) 

23.1(3.1) 

23.9(1.4) 

25.1(2.4) 
24.7(3.4) 
22.8(4.5) 



Table 2 

Results for the conversion factor F for various 
scales Q' and fixed Q — r^ 1 at zero lattice spac- 
ing. Also shown are the estimates for mi^ ex- 
pressed through the common scale r^ 1 in MeV 
(MS-scheme at /j, = 2 GeV). 



Ref. 


(Q') 
nis 


Q' 


F 


mi ro) 




13 


106(7) 


m p 


0.90(1) 


95(6) 




n| 


114(2)(1) 


m p 


0.86(2) 


98(2)(I) 




13 


106(2)(8) 


niK* 


0.87(3) 


92(2)(7) 




T§ 


97(4) 


F K 


1.02(2) 


99(4) 



to determine the ratio (Q/Q')i a t in the continuum 
limit for Q = r^ 1 and Q' = m p , m,js*,Fjc. The 
factor F can then be computed using the phe- 
nomenological values of ro = 0.5 fm and Q' . The 
results are shown in Table @. The first observa- 
tion is that F varies by up to 15% for the var- 
ious scales that are compared. Once the results 
have been converted to the common scale ro , the 
estimates for m s in the continuum limit show re- 
markable consistency. This is actually surprising, 
since the various simulations differ significantly, 
not only in terms of the fermionic discretisation, 
but also by the chosen method to relate the lat- 
tice estimate to the MS scheme (mean-field per- 
turbative, non-perturbative via RI/MOM or SF). 
Lastly, I want to stress that a meaningful com- 
parison of this kind can only be made in the con- 
tinuum limit. 



2.2. Ginsparg- Wilson fermions 

Lattice actions that satisfy the Ginsparg- 
Wilson relation ||l| and thus preserve chiral sym- 
metry at non-zero lattice spacing are now rou- 
tinely used to compute many phenomenologically 
interesting quantities. The most widely used im- 
plementations are based on overlap ||| or do- 
main wall fermions |23| . More recently, results 
for fixed-point (FP) J24j or chirally improved (CI) 
|25|] actions, which both provide approximate so- 
lutions to the Ginsparg- Wilson relation, have also 
become available (see [Q for a review). 

The mass of the strange quark was actually one 
of the first quantities computed using domain wall 
fermions |27| . Since then many systematic effects 
have been studied: discretisation errors have been 
estimated by comparing results at different lat- 
tice spacings, although no continuum extrapola- 
tions have been performed so far. Quenching ef- 
fects have not been investigated, since the numer- 
ical effort to simulate quenched Ginsparg- Wilson 
fermions is already comparable to a dynamical 
simulation using Wilson quarks |p8| . 

We have seen that non-perturbative renormal- 
isation is an important ingredient in order to en- 
hance the credibility of lattice estimates for quark 
masses. The RI/MOM prescription has already 
been applied for both domain wall [^9| and over- 
lap fermions [^0j. The Schrodingcr functional, 
due to its inhomogeneous boundary conditions is 
somewhat harder to realise for Ginsparg- Wilson 
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fermions (an attempt has been made in ref. [pl| in 
the case of domain wall fermions). The authors of 
J32J have proposed a general strategy to compute 
the renormalisation factors of quark bilinears for 
overlap fermions in the SF via an intermediate 
Wilson regularisation. This procedure avoids the 
direct formulation of the SF for Ginsparg- Wilson 
fermions, at the expense of sacrificing the abil- 
ity to predict the quantity used to match to the 
Wilson results in the intermediate step. 

Results for m s and m s /fh obtained using ei- 
ther domain wall or overlap fermions are com- 
piled in Table |^. The numbers for m s are broadly 
consistent among each other, and also with the 
continuum estimates discussed earlier. However, 
the overall uncertainties are still somewhat larger 
compared with Wilson or staggered fermions. In 
most cases, the quoted systematic uncertainties 
are dominated by lattice artefacts, which have not 
been removed by performing continuum extrap- 
olations. However, the overall dependence of m s 
on the lattice spacing appears to be fairly weak. 
These indications of good scaling properties of 
Ginsparg- Wilson fermions must be corroborated 
with more statistics and systematic investigations 
of renormalisation and finite-a effects. 

2.3. Dynamical simulations 

The most important challenge for current lat- 
tice studies is surely the quantification of quench- 
ing effects. Perhaps the most comprehensive 
study so far has been published by CP-PACS 
some time ago [[[3). They investigated Nt = 2 
flavours of dynamical quarks, using mean-field 
improved Wilson fermions and the Iwasaki gauge 
action at three different values of the lattice spac- 
ing in the range a rj 0.1 — 0.2 fm. Input val- 
ues for the dynamical quark masses were chosen 
such that TOps/my = 0.58 — 0.8. Estimates of 
m s based either on vector or axial vector Ward 
identities (VWI and AWI respectively) can be ex- 
trapolated to a common value in the continuum 
limit. By comparing their results to quenched 
data, CP-PACS find that dynamical simulations 
lead to a decrease in the value of m^ s by about 
25%. The dynamical data also appear to be less 
sensitive to the choice of quantity used to fix the 
bare strange quark mass (ttjk or m$). The main 



findings of CP-PACS are summarised in Fig. [|, 
as well as in the estimates 



mf b (2 GeV) = MeV, m s /m w 26. 



(7) 
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Figure 1. Scaling plot of dynamical and quenched 
estimates for m s by CP-PACS @. 



While CP-PACS's results represent a major 
step forward, the limitations of current dynamical 
simulations are apparent. In particular, smaller 
dynamical quark masses should be simulated in 
order to quantify quenching effects more reliably. 
Moreover, simulating two flavours of dynamical 
quarks still does not represent the physical situ- 
ation. Discretisation errors must be investigated 
more closely, by simulating smaller lattice spac- 
ing and comparisons with alternative fermionic 
discretisations. Finally, non-perturbative renor- 
malisation should be implemented. Attempts to 
improve the situation in each of these areas have 
been reported at this conference. 

The JLQCD and UKQCD Collaborations have 
both used non-perturbatively 0(a)-improved 
Wilson fermions for N t = 2 @ at j3 = 5.2. 
Whereas JLQCD uses sea quark masses cor- 
responding to tops/tov = 0.6 — 0.8, UKQCD 
simulates even lighter masses, thereby running 
the risk of suffering from finite-volume effects. 
JLQCD have reported finite-size effects at the 
3 — 5% level only at their lightest quark mass. 
They also extrapolate their results in m sc ' A to the 
physical u and d-quark masses. The quark masses 
are defined through the PCAC relation (AWI), 
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Table 3 

Results for m s and the ratio m s /m, obtained using domain wall (DW) and overlap (Ov) quarks. Con- 
ventions are identical to Table 0. 



Collaboration 


Rcf. 


action 




NP-Rcn. 


m s [MeV] 


Ul s f TTL 


RBC 


( 




1 


DW/DBW2 
DW/W 


-0.099 
-0.099 


RI 


133(3) 
126(3) 


~ 26 


CP-PACS 


34 


1 


DW 


>0.066 




99(2)(6) 


26.3(2.3) 


RBC 


1 


i 


1 


DW 


0.093 
0.123 


RI 


105(6)(21) 
100(5) (20) 




BSW 


I 




1 


DW 


>0.073 




96(26) 




DcGrand 


1 


§ 


1 


Ov 


0.090 
0.125 




110(7) 
105(7) 


24.41(5) 
24.40(4) 


C+H 


37 


1 


Ov 


0.147 




115(8) 


21.9(2) 


GHR 


38 




Ov 


0.086 


RI 


102(6)(18) 




HJLW 






1 


Ov 


0.093 
0.123 


SF 


96(5)(7) 
99(4)(7) 


26.1(3.0) 
25.9(3.4) 



and they expect some systematic effects to can- 
cel, since the AWI mass is defined through a ratio 
of matrix elements. The results for m s and m are 
20 — 30% smaller than in the quenched approx- 
imation. Scaling violations appear to be small, 
as there is good agreement with the previous CP- 
PACS results Jl3| , but a continuum extrapolation 
is still lacking. 

UKQCD perform extrapolations in the valence 
quark mass at fixed values of m soa and then mon- 
itor the results for decreasing sea quark mass, in 
order to search for trends which indicate the pres- 
ence of dynamical quark effects. Unlike JLQCD 
they define the quark mass via the vector Ward 
identity. UKQCD's results show no clear signs of 
dynamical quark effects: although there is some 
trend in the data as m sca is decreased, it is sta- 
tistically not significant, and the same is true for 
the comparison with quenched data at a matched 
value of the lattice spacing. Their result for the 
ratio m s /rh is consistent with ChPT. 

Hein et al. jl3) have used dynamical configu- 
rations provided by the MILC Collaboration Jl4| 
for 2+1 flavours of improved staggered quarks. In 
this formulation the large flavour-changing inter- 
actions, which are typically encountered for con- 
ventional staggered quarks are suppressed due to 
the use of fat links. The lightest quark mass in 
(44| corresponds to mps/mv ~ 0.37, at a lattice 
spacing of a ~ 0.13 fm. Mean-field improved per- 



turbation theory was used when matching to the 
MS-scheme; the one-loop coefficients were found 
to be of the same order of magnitude as for Wil- 
son fermions. Results for m s + m obtained for 
Nf — 0, 2 and 2+1 dynamical flavours show that 
the large dependence on the quantity that sets 
the lattice scale is greatly reduced in the dynam- 
ical case. Low values for the light quark masses 
are typically preferred; for Nf = 2 + 1 flavours 
Hein et al. quote 

(m s +m) m (2 GeV) = 78 + 14 MeV, (8) 

where the scale has been taken from the K* — K 
mass splitting. The quoted error includes the un- 
certainty from neglecting the two-loop term in 
the perturbative matching procedure, which is 
estimated to be as large as 20%. The value of 
(m s + fh) in the dynamical case is lower by 15 - 
20% compared with the corresponding quenched 
result. For a better understanding of the system- 
atics in the improved staggered formulation, it 
would, in my opinion, be extremely helpful if a 
quenched value for m s were available in the con- 
tinuum limit. 

The ALPHA Collaboration has reported re- 
sults on the non-perturbative running of the 
quark mass in the SF scheme for the 2-flavour 
case Q. This is one part of the two-step pro- 
cedure, which will ultimately relate the bare cur- 
rent quark mass on the lattice to the renormali- 



G 



sation group invariant (RGI) quark mass M 
In Fig. |^ the running mass in the SF scheme 
is plotted as a function of the renormalisation 
scale. One sees that the perturbative evolution 
follows the numerical data down to fairly small 
values of the scale. At its lowest value, i.e. at 
ji = l/(2L max ), ALPHA obtain their preliminary 
result for the matching of the running mass in the 
SF scheme to the RGI quark mass: 

M 

— = 1.236(15), ln(AL max ) = —1.85(13). (9) 
m S F 

The scale L max must still be related to some phys- 
ical quantity. Further studies will include data at 
larger values of L/a, so that lattice artefacts can 
be eliminated completely. 



m( M )/M . 




Figure 2. Running quark mass in the SF scheme 
for N f = 2 §|. 

2.4. The charm quark mass 

The charm quark is too heavy to be described 
by ChPT, and too light for an efficient non- 
relativistic treatment. In two recent projects 
JI(i| , [47t the charm quark mass was computed us- 
ing the relativistic formulation. In this case, the 
issue of controlling lattice artefacts is even more 
important as in the light quark sector. The AL- 
PHA Collaboration |l7j have used 0(a) improved 
Wilson fermions in the quenched approximation 
at four (3- values. They employ several different 
definitions of the RGI charm quark mass, e.g. 

r M c = Z M r m c [1 + (b A - b P )am qtC ] (10) 
r M c = Z M Zr m q:C [1 + b m am qtC ] , (11) 



where m c and m q:C denote the bare masses de- 
fined through the axial and vector Ward iden- 
tities, respectively. Non-perturbative values of 
renormalisation factors and improvement coef- 
ficients as reported in [{spsfl have been used 
throughout, in order to guarantee a controlled 
extrapolation to the continuum limit. The input 
quantities which fix the bare charm quark mass in 
terms of meson masses were mD B , tuk and the ra- 
tio m s /fh from ChPT. In this way chiral extrapo- 
lations can be avoided completely. By demanding 
that different definitions of the RGI charm quark 
mass extrapolate to a common continuum limit, 
they find r M c = 4.19(11) at zero lattice spacing, 
which translates into 

m^(m c ) = 1.301(34) GeV, (12) 

where the 4-loop RG functions in the MS scheme 
have been used. This result agrees well with the 
estimate by the SPQ c dR Collaboration p6{ 

rrjfS(m c ) = 1.26(4)(12) GeV, (13) 

which has been obtained at a 0.07 fm. 
2.5. Summary: Quark masses 

A detailed comparison of quenched data shows 
that the mass of the strange quark is probably 
the most precisely determined quantity in the 
quenched approximation: for a given choice of 
lattice scale, the relative accuracy in the contin- 
uum limit is about 5%. If the if -meson is used to 
fix the bare ra s , the final result varies between 95 
and 115 MeV, depending on the choice of lattice 
scale. There is a 20% increase in m s if the mass 
is fixed through a vector meson like the (j>. 

It is now clear that current lattice estimates 
for the strange quark mass are completely dom- 
inated by quenching effects. First attempts to 
quantify the quenching error indicate a decrease 
in m s of about 20%, but simulations at smaller 
quark masses and lattice spacings are required to 
confirm this. 

Estimates for the ratio m B /m based on chiral 
extrapolations tend to agree with ChPT, both for 
quenched and unquenched simulations. 

Finally, non-perturbative renormalisation and 
O(a) improvement enable one to perform con- 
trolled continuum extrapolations of the charm 
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quark mass, thereby yielding precise quenched es- 
timates. 

3. QUENCHED CHIRAL LOGARITHMS 

In this section we will be concerned with 
quenching artefacts in the chiral expansion. It 
is well known that the quark mass behaviour of 
hadron masses and matrix elements is modified 
in the quenched approximation The 
main qualitative difference is that flavour sin- 
glets do not decouple in the quenched version of 
the effective low-energy theory. This is because 
most of the loop contributions which lead to the 
decoupling of flavour-singlet fields are removed, 
and hence the latter have to be treated on the 
same footing as the octet fields. Incorporating 
the flavour-singlet explicitly into the chiral La- 
grangian introduces new LECs: one is the mass 
scale m of the flavour singlet, the other is a$, 
which multiplies the kinetic term in the flavour- 
singlet part of the effective Lagrangian: 



(14) 



The flavour-singlet propagator develops a double 
pole in the quenched theory, which is propor- 



tional to 



Since mo does not vanish in the 



chiral limit, this gives rise to a new type of chi- 
ral logarithm, which diverges as the octet mass 
vanishes. The modified chiral expansion for the 
pseudoscalar octet mass at NLO thus reads 



= B {l-(*-!^)(ln y + l) 
+ [(2a 8 - as) - ^a®] y}. 



(15) 



For notational convenience we have introduced 



(47rF ) 2 7V f ' V (4ttF ) 



2B m 



(16) 



and rescaled the LECs Li to a; = 8(4-7r) 2 Li. A 
phenomcnological estimate of 5 is obtained from 
the relation 



mr, 



TO 2 , + TO 2 — 2TO-K 



(17) 



which gives 5 w 0.18. There are several meth- 
ods which allow the determination of 8 and a$. 
The most straightforward is to study the quark 
mass dependence of TOp S and compare it with 



eq. (y_5J). Another method is to compute the ra- 
tio of flavour-singlet and flavour-octet contribu- 
tions to the quenched rj' correlator. The former 
are the so-called "hairpin diagrams" . Finally, the 
Witten-Veneziano formula |52| provides a link be- 
tween the flavour-singlet mass too and the topo- 
logical susceptibility \t- 



(18) 



The effects of quenched chiral logarithms are ex- 
pected to show up very near the chiral limit. 
However, in this regime quenched QCD is af- 
flicted with the occurrence of so-called "excep- 
tional configurations" , i.e. unphysical near-zero 
modes, if conventional discretisations such as Wil- 
son fermions are used. Recently this has stimu- 
lated the use of discretisations which do not suffer 
from this problem, such as overlap and domain 
wall fermions, fixed point actions and twisted 
mass QCD (tmQCD) |§. Also the "Modified 
Quenched Approximation" |54|j in conjunction 
with Wilson fermions has been applied. In addi- 
tion, one has to avoid finite-volume effects when 
going very near the chiral limit, which may fake 
the signature of the quenched chiral logs. It is 
also clear that there is a special role for discreti- 
sations which preserve chiral symmetry, for which 
the Witten-Veneziano formula is exact J55| , p6[ . 
Moreover, the comparison of lattice data with the 
predictions of quenched lattice QCD at non-zero 
lattice spacing is justified in this case. 

Lattice data for to ps show very low sensitivity 
on a$ in general. The RBC Collaboration finds 
]5S| ] that in the studied mass range the combina- 
tion a^TOpglnTOpg (or a$>y\ny in eq. ([D])) stays 
practically constant, so that there is hardly any 
sensitivity to the value of a$. CP-PACS report 
|l^1 that fits in which a$ is treated as a free pa- 
rameter are either unstable or yield results that 
are consistent with zero. Thus, in most studies it 
is assumed that a$ = 0. 

Table [| shows a compilation of recent determi- 
nations of S. Overall one observes large variations 
in the current estimates for <5, with results from 
simulations employing conventional fcrmionic dis- 
cretisations being somewhat lower that those us- 
ing Ginsparg- Wilson fermions (exceptions are the 
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Table 4 

Recent determinations of S. The second and third columns show the minimum (bare) quark mass in McV 
and the smallest pion mass in units of the spatial lattice size, respectively 



Collaboration 


Rcf. 


m q [MeV] 


m™ m L 


a [fm] 


action 


5 


CP-PACS 


I 


12 


1 


20 


6 


0.05-0.1 


W 


0.10(2) 


QCDSF 


I 


57 




22 


4.7 


0.05-0.1 


sw 


0.14(2) 


FNAL 


I 


54 




14 


3 


0.17 


sw 


0.065(13) 


MILC 


I 


14 




15 


3.4 


0.13 


KS 


0.061(3) 


Kentucky 


I 


59 




19 


3.2 


0.157 


Ov 


0.2-0.4 


RBC 


I 


58 




28 


3.8 


0.104 


DW 


0.05(2) 


DeGrand /Heller 


I 


5G 




30 




0.13 


Ov 


0.093(28) 


Chiu+Hsieh 


I 


60 




80 


2.5 


0.147 


Ov 


0.203(14) 


Bern 




Cil 




38 


2.6 


0.13 


FP 


0.23(7) -0.30(18) 


Kentucky 


I 


32 




14 


2.9 


0.20 


Ov 


0.26(3) 


BGR 


1 


i 


1 




2.7 


0.16 


FP 


0.17(2) 










2.9 


0.15 


CI 


0.18(3) 



low values quoted in refs. |5^.|5(|). Despite the 
fact that the sensitivity to quenched chiral logs 
is still quite low, there is mounting evidence for 
a non-zero value for S. In fact, more recent sim- 
ulations, which employ lattice chiral symmetry, 
produce results that are compatible with the phe- 
nomenological estimate based on eq. ([l7|). How- 
ever, a remaining matter of concern is the pos- 
sibility that finite-volume effects may distort the 
mass dependence towards the chiral regime. This 
must be addressed in simulations on larger vol- 
umes. So far no convincing signal for a$ has been 
observed, mainly because the sensitivity of lattice 
data is even worse than for 6. 

4. IS THE UP-QUARK MASSLESS? 

Let us now return to the Kaplan-Manohar am- 
biguity. The fact that the effective chiral La- 
grangian beyond leading order is invariant under 
simultaneous transformations of the mass matrix 
and a subset of LECs implies an uncertainty in 
the size of the NLO correction to the mass ratio 
m u /md. This correction, Am, is given by 



A M = 



(m| 



ml) 



(2a$ — as) + chiral logs. (19) 



While as can be estimated from the ratio Fk/F^, 
there exists no experimental information which 
would allow an unambiguous determination of ag 
or indeed the linear combination 2ag — |64|J65| . 



The value of as and thus Am must then be fixed 
by invoking (plausible) assumptions beyond chi- 
ral symmetry considerations. Proceeding in this 
way Leutwyler [M obtained 



< A M < 0.13. 



(20) 



This result excludes the possibility that m u /md = 
0, which is only possible for a large negative 
value for Am- Therefore, a massless up-quark, 
which presents a simple and elegant solution to 
the strong CP problem is an unlikely scenario. 
However, an estimate of Am based on first prin- 
ciples is still lacking. Given the importance of 
the strong CP problem, it is desirable to tackle 
this question by means of lattice determinations 
of 2a8 — as- 

Here, partially quenched simulations of lattice 
QCD play an important role: for unequal sea 
and valence quark masses the effective chiral La- 
grangian is parameterised in terms of the same 
LECs as the physical theory. Thus, as long as the 
correct number of dynamical quark flavours (i.e. 
Nf — 3) is used, simulations based on unphysi- 
cal mass combinations provide phenomenological 
information |66|| , provided that the regime of va- 
lidity of ChPT and the range of simulated masses 
overlap. 

A general method which allows for the extrac- 
tion of LECs with good statistical accuracy has 
been introduced by ALPHA [p7| . Since the ex- 
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pressions of ChPT are valid for arbitrary quark 
mass, one can introduce a reference value m re f 
and define the ratios 

R F {x)=F PS (m)/F PS (m ref ) (21) 



Rm(x) = 




where x is m/m re f. The ratios Rf, Rm are uni- 
versal functions of the dimensionless mass param- 
eter x. All renormalisation factors drop out, so 
that these quantities can be extrapolated to the 
continuum limit straightforwardly. Due to the 
correlations between the numerator and denom- 
inator one can expect good statistical precision 
in numerical evaluations. The LECs at NLO can 
then be determined by comparing numerical data 
for i?p , Rm with the corresponding expressions in 
quenched H@|l) or partially quenched p8| , j39[ 
ChPT. 




- 5 = 0.12, a t =0.0, (2a 8 -a 6 )<°)=0. 345(g)- 




0.5 1 1.5 



X 

Figure 3. Results for Rf and Rm in the quenched 
approximation |57j. 

ALPHA have tested their method in the 
quenched approximation for m re f « m s , and x in 
the range 0.75 < x < 1.4 ^j?]]. Their results are 
shown in Fig. [|. The ratio Rp, which is predicted 
by quenched ChPT to rise linearly with the quark 
mass, is modelled very well by the data, result- 
ing in stable estimates for a§. Numerical data 
for the ratio Rm, however, show an almost con- 
stant behaviour, while quenched ChPT predicts 



the presence of linear terms, as well as chiral logs. 
Thus, if the range of simulated masses lies indeed 
within the regime where ChPT at NLO is valid, 
then the constant behaviour of the data must be 
the result of a strong cancellation between the 
various contributions at NLO. 

Bardeen et al. |34| have performed global fits 
to pseudoscalar correlation functions to extract 
masses and various LECs. Using the pole shifting 
prescription of the "Modified Quenched Approx- 
imation" (MQA) at (3 — 5.7 and mean-field im- 
proved Wilson fermions, they report an estimate 
for 05 which is about 3 times larger than AL- 
PHA'S. A large fraction of this discrepancy can 
be attributed to the fact that renormalisation fac- 
tors for the axial current have not been taken into 
account |7^]. In fact, if the data for the decay con- 
stant published in |34j are subjected to ALPHA'S 
ratio method, then the estimates for a§ are con- 
sistent among the two studies. Thus, the source 
for the original discrepancy is a combination of 
renormalisation effects and lattice artefacts, with 
the ratio method being a lot more stable against 
these systematic effects. 

The ratio method has also been applied in two 
recent simulations of partially quenched QCD. 
UKQCD jnj have simulated two flavours of 0(a) 
improved dynamical Wilson fermions at a = 
0.1 fm and a fixed value of the sea quark mass of 
TO sea ^ o.7m s , which corresponds to mps/my ~ 
0.58 flr2| . The mass dependence of the ratios Rf 
and Rm was subsequently studied as a function 
of the valence quark mass at fixed m soa . The 
overall behaviour of Rf and Rm was found to be 
similar to what is observed in the quenched case. 
However, the strict linearity of Rf for Nf = is 
modified by a slight curvature, which may signal 
the expected presence of chiral logs. UKQCD's 
results for the LECs are listed in Table ||. 

The Ohio group [|7J| have reported results for 
LECs from simulations using N{ = 3 flavours 
of dynamical staggered quarks. The sea quarks 
are somewhat lighter than in UKQCD's study, 
while the lattice spacings are larger (a rs 0.15 
and 0.28 fm). The dominant systematic effects 
are large flavour-symmetry violations, which are 
estimated by applying hypercubic blocking to the 
gauge configurations before the computation of 



10 



Table 5 

Results for the LECs at NLO for several values of Af. The last two lines display the results for the 
"standard" phenomenological estimates, as well as the numbers required for the up-quark to be massless. 



Collab. 


Rcf. 


Nt 


«5 


a$ 


2a 8 - a 5 


comments 


ALPHA 


§§ 




"3" 


0.99(6)(20) 
0.75(6)(20) 


0.50(4)(20) 
0.45(4)(20) 


0.35(5)(15) 
0.15(5)(15) 


S = 0.12, a$ = 


UKQCD 





2 

"3" 


1.22(13)(25) 
0.98(13)(24) 


0.79(6)(21) 
0.59(6)(21) 


0.36(10)(22) 
0.20(11)(22) 




osu 


n 


3 






0.33(2)(17) 








3 
3 


0.5(6) 
0.5(6) 


0.76(4) 
-0.9(4) 




"standard" 
m u = 



observables. The value of 0:5 differs significantly 
when it is evaluated on blocked and unblocked 
configurations. In the future this issue will be ad- 
dressed by comparing Monte Carlo data to ChPT 
for which the effects of flavour-symmetry break- 
ing are included Q . 

From the compilation of results in Table || one 
concludes that the Af-dependence of LECs is 
quite weak. If one uses the expressions of ChPT 
for Af = 3, even though actually N{ = or 2 was 
used to generate the data, one still observes con- 
sistency within errors (cases labelled Af ="3" in 
Table ||). A comparison with the "standard" val- 
ues for the LECs estimated in the continuum, as 
well as with those values that would be required 
to support the notion of a massless up-quark, 
shows that the latter scenario is strongly dis- 
favoured by lattice calculations: the quark mass 
behaviour of Rm would have to be radically dif- 
ferent to accommodate a large negative value of 
ag- In short, calculations based on first principles 
support the analysis of ref. (20| . 

The most important issue at this stage is 
whether or not the quarks used in the simula- 
tions arc light enough to justify the comparison 
with ChPT. This has been addressed in several 
contributions to the panel discussion on chiral ex- 
trapolations at this conference Q. Although lat- 
tice estimates for some of the LECs make sense 
phenomenologically (e.g. UKQCD's result for a$ 
in fn^ is consistent with the experimental value 
of Fk/Fk J7l[]), there is evidence that the depen- 
dence of m,p S and Fpg on the sea quark mass is 



not in agreement with ChPT if the sea quark 
mass corresponds to ?7ips = 550 — 1000 MeV 
{TaJI. This affects attempts to perform extrap- 
olations in m sea to the physical point defined by 
fnps/ m v ~ 0.17. A similar study j77) concluded 
that sea quark masses of current simulations and 
ChPT at NLO overlap only marginally. 

In addition to many numerical results, there are 
also new analytic developments. Aubin et al. |74| 
have analysed the structure of chiral logarithms 
under staggered flavour-symmetry breaking for 
Af = 2 + 1 flavours. The basic idea is to generalise 
the Lagrangian of Lee and Sharpe j78| , which de- 
scribes a single staggered field up to 0(a 2 ), to the 
3-flavour case. Normally, each staggered flavour 
describes four internal fermions (called "tastes" 
in {74]]). In order to obtain one taste per flavour 
one then takes the fourth root of the determi- 
nant. For the actual calculation of chiral logs in 
the 2+1 flavour case, this implies that one starts 
with 4+4 "tastes" and then cancels the unwanted 
loop contributions by multiplying them with sim- 
ple weight factors. The resulting expressions for 
pseudoscalar meson masses fit the data by MILC 
p4| very well, unlike the corresponding formulae 
in "ordinary" partially quenched ChPT j68pS[ |. 

In a recent paper J79| Rupak and Shoresh have 
incorporated the explicit chiral symmetry break- 
ing of Wilson fermions directly into the effective 
Lagrangian. The starting point is Symanzik's ef- 
fective action. Including all operators up to di- 
mension 5, they find that chiral symmetry is bro- 
ken not only by the quark density i^ip, but also 
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by the Pauli term ^(t^F^i/j, which is multiplied 
by the coefficient c sw . Rupak and Shoresh then 
consider simultaneous chiral expansions in the pa- 
rameters 

2B a m 2W ac sw 



A 2 



A 2 

x 



(23) 



The Pauli term generates additional operators 
in the effective chiral Lagrangian, starting at 
leading order. Consequently, there are a num- 
ber of additional LECs. Up to NLO these are 
Wo, W4, . . . , Ws, where the subscripts are chosen 
in analogy with the operators appearing in the 
conventional Lagrangian. The resulting expres- 
sions for pseudoscalar masses and decay constants 
in the partially quenched case may be used to ex- 
tract the LECs at small but non-zero lattice spac- 
ing. 

5. SUMMARY 

Lattice determinations of quark masses have 
entered a mature stage, where the dominant sys- 
tematic error is quenching. The quantification of 
dynamical quark effects will be the main subject 
of investigation in the coming years. The com- 
parison of lattice data with effective low-energy 
theories such as ChPT has turned into a major 
activity. Here, the goal is not just to verify ChPT 
but to exploit its predictions in order to make con- 
tact with the regime of very light quarks, which is 
difficult to access directly in simulations. So far it 
is not entirely clear whether the sea quark masses 
used in simulations are small enough to justify ex- 
trapolations to the chiral regime based on ChPT 
at NLO. Simulations of partially quenched QCD 
show that the scenario of a massless up-quark is 
not reflected at all in the valence quark mass be- 
haviour of pseudoscalar quantities. Further simu- 
lations with smaller quark masses should be per- 
formed to corroborate these findings and to set- 
tle the issue of the mass range in which ChPT at 
NLO can be expected to be valid. 
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